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ABSTRACT 


When a circular cylinder is inmersed in a fluid of infinite extent 
and oscillated in such a manner that drag effects are negligible, it is 
well known that the closed form solution gives an added mass coefficient 
of 1.0. However, when boundaries are placed in proximity to the cylinder, 
and the cylinder oscillated parallel to the boundary, the added mass 
varies from the classical result. The forces acting on the cylinder in 
the latter case are dependent upon the proximity of the boundaries as 
well as the frequency of oscillation. 

It was the purpose of this study to experimentally determine the 
total force on such a cylinder in proximity to an impermeable bottom 
surface and under a free surface and to relate this force to the added 
mass coefficient. Experimental data was reduced using a computer program 
and the results presented graphically as a function of dimensionless 


parameters. 
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NOMENCLATURE 
Description 

Radius of Circular Cylinder 
Relative Depth of Sukmergence 
Depth of Submergence of Cylinder 
Total Maximum Force Per Unit Length 
Acceleration of Gravity 
Relative Water Depth 
Water Depth 
Phase Angle 
Waveheight 
Relative Displacement 
Half Amplitude of Cylinder Oscillation 
Dimensionless Frequency Parameter 
Water Density 
Circular Frequency of Cylinder Oscillation 


Water Viscosity 





DIMENSIONLESS NOTATION 


F/pgra%o 2X Force Coefficient 

ofa/g =v Frequency Parameter 

aya ="d Relative Depth of Submergence 
ma = h Relative Water Depth 

e/a =9x Relative Displacement 


2paXs/\ Reynold's Number 
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I. INTRODUCTION 


During the past two decades, many new submarine sitructures, such as 
pipelines, oil storage tanks and large aquaducts have been either sunk 
or moored on the continental shelf. Consequently, the calculation of the 
forces exerted on such objects is of practical importance. 

In an early paper [1], Stokes showed that the expression for the 
force on a pendulum-like system oscillating in an infinite fluid was both 
acceleration and velocity dependent. 

In more recent years, Morison, et al, [2], from an investigation of 
wave forces exerted on piles, showed that when the characteristic dimen- 
Sion of a body, i.e., pile diameter, was small compared to wave length, 
the force consisted of both a drag and an inertial teym. However, the 
coefficient of each of these terms depends upon the geometry of the 
object in question, the amplitude of the fluid motion, the proximity of 
a free surface, and the proximity of an impermeable bottom surface. 

Several investigations have been conducted since Morison's work, which 
have dealt with sinusoidal flow about various objects. Keulegan and 
Carpenter [3] performed an experimental study with cylinders and plates 
in a rectangular tank in which a standing wave was generated in water. 
Forces on the submerged bodies were measured with a dynamometer. Inertial 
and drag er icin were correlated with a dimensionless "period 
parameter" for both unseparated and separated flow, the period parameter 


being proportional to the relative displacement of the surrounding fluid. 
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Heinzer [4] carried out a visual study of the flow field character- 
istics of a cylinder sinusoidally oscillated in a still fluid. Flow 
patterns included laminar boundary layer, symmetric separation, and 
Karman vortex street. Heinzer's major conclusion was that wake charac- 
teristics were primarily dependent upon the relative amplitude of oscil- 
lates. A consideration of this type was used in the study which follows. 
Small amplitude oscillations were utilized at various frequencies to 
insure unseparated flow. 

Sarpkaya and Garrison [5] studied constantly accelerated, 
unidirectional flow past a circular cylinder. Forces were measured 
experimentally and the nature of the flow determined visually. Force 
coefficients were calculated as functions of the relative displacement 
of the cylinder. It was shown that the total resistance as well as the 
drag and inertial forces are representable as functions of the relative 
displacement of the fluid. 

Ward and Dalton [6] analyzed the forces exerted on a cylinder ina 
purely sinusoidal flow field in which synmetrically located vortex pairs 
were developed. Forces acting on the cylinder were calculated using 
Bryson's [7] potential flow model including a pair of vortices but no 
attempt was made to verify these forces experimentally. 

Later Dalton and Hamann [8] conducted an experimental investigation 
of a sinusoidally osciallating cylinder in a quiescent fluid in an attempt 
to represent a Sinusoidal wave moving past a single, fixed platform leg. 


The conclusions drawn from their investigation were mostly qualitative. 
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The purpose of this study was to provide a correlation between 
the added mass coefficient and appropriate dimensionless parameters. A 
right circular cylinder was oscillated in a quiescent fluid at ampli- 
tudes which insured unseparated flow. The results are presented in 
dimensionless form as a function of a dimensionless frequency parameter, 
relative water depth, relative depth of submergence, and relative 


amplitude of oscillation. 
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IT. THEORETICAL CONSIDERATIONS 


A schematic representation of the problem is shown in Figure (1). A 
right circular cylinder of radius a was submerged to a depth d in water 
of depth h. The cylinder was oscillated sinusoidally at a circular fre- 
quency o with half amplitude of oscillation x Of primary interest was 
the total maximm force per unit length F exerted by the water on the 
cylinder, the phase of this force and the amplitude of the generated waves. 

The camlete theoretical solution to the problem is complicated, even 
when the amplitude of oscillation is kept small enough to justify an 
assumption that the flow remained unseparated. Therefore, in order to 
make sane headway a dimensional analysis was applied to the problem. It 
is known "a priori" that the maximum force per unit length of cylinder has 


the following functional dependence. 


F = fy (invelyai,c X72 Ae raly) (1) 


where h denotes the water depth, d the depth of submergence, a the cylinder 
radius, o the circular frequency of the imposed motion, X the half ampli- 
tude of the imposed motion, p the water density, g the acceleration of 
gravity, and » the water viscosity. 

Dimensional analysis of the physical quantities of Equation (1) 
resulted in the following dimensionless groups 
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where = (1S the well known Reynold's number representing the ratio of 
inertial to viscuous forces. Since the cylinder was oscillated at small 
amplitudes it was assumed that flow around the cylinder was unseparated 
and viscous forces small. Therefore, the Reynold‘s number was disregarded 


and Equation (2) written as 


2 
F/ogna’s x = £0, : (3) 


| [Oy 
m1] 
m {1X 


A second reason for ignoring a dependence upon Reynold's number was 
provided by Heinzer [4]. His visual studies showed that wake characteris- 
tics were primarily dependent upon the relative amplitude of oscillation 
and that Reynold's number effects were not evident. 

As long as the relative displacement is small, the flow remains 
essentially unseparated and, accordingly, the force is also independent 
of the relative amplitude. The independence of the force from this para- 
meter for small values was demonstrated experimentally by Keulegan and 
Carpenter [3]. In the present ee cero no attempt was made to study 
the effects of this parameter; consequently, it was held fixed at - Oral 
throughout the experiment. 

The objective of this experimental work was to show a correlation 
between the dimensionless force coefficient on the left hand side of 


Equation (3) and the dimensionless parameters on the right hand side. 
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IIIT. EXPERIMENTAL INVESTIGATION 


A. APPARATUS 
1. Wave Channel 

The wave channel was designed and built by Shiller [9] for an 
earlier thesis. The channel was designed with a rectangular cross section 
two feet deep and fifteen inches wide. It was fabricated in sections, 
each section consisting of two vertical walls and one horizontal floor 
member. The vertical walls were stiffened by use of three four foot long 
2 x 4's attached to 3/4 inch x 8 foot x 3 foot sheets of plywood at the 
ends and middle of the longest dimension. The channel floor member was 
composed of two 15 inch x 8 foot long plywood sheets separated by a 
1 inch x 4 inch strip around the edge. The bottom and side members for- 
formed the basic channel cross section. The bottom and side members were 
mated together with 3/8 inch threaded rods. Six of these 8 foot sections 
were sealed and bolted together to forma 48 foot long channel. Square- 
ness of the channel cross sections was maintained by use of spacers placed 
between the stiffened members near the floor and across the top. 

One end of the channel was fitted with a wave generating mechanism 
which was not used in this investigation. 

At a location 30 feet from the wave generator end of the channel, 
a three foot wide plexiglass observation window extending the full height 
of the channel was installed in both walls as shown in Figure (2). The 
test cylinder was sinusoidally oscillated at a position near one end of 


the plexiglass window. 
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The generated waves were dissipated on a beach structure located 
at each end of the channel. The beach as shown in Figure (3) consisted 
of three layers of Peeeor ata stainless steal sheeting. Layer separation 
and structural strength were provided by 2 x 4 runners with their largest 
cross sectional dimension perpendicular to the sheeting. The slope of 
the beach was 6:1. 

2. Carriage and Driving Mechanism 

The driving mechanism as shown in Figure (4) consisted of a vari- 
able speed transmission driven by a two-horsepower inductance motor. The 
motor and transmission were mounted on a 3/4 inch steel plate which was 
positioned and bolted to the top of the channel. The plate was further 
braced using 1/2 inch angle iron supports which were bolted across the 
top of the channel and diagonally to the channel side wall. The variable 
Speed transmission provided an output shaft speed range extending fram 22 
to 180 rem which in turn provided dimensionless frequency parameters of 
fram about 0.04 with a six inch diameter cylinder to about 1.7 witha 
four inch diameter cylinder. 

The output shaft of the transmission was fitted with a face plate. 
The face plate as shown in Figure (5) was connected to the drive rod 
through a bolt and slot assembly which provided a drive rod stroke which 
could be varied between zero and six inches. 

The carriage as shown in Figure (6) provided a movable mounting 
platform capable of being positioned over the top of the wave channel. 
The carriage was driven sinusoidally by the driving mechanism at 


amplitudes of about 10% of the cylinder diameter. 
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PIGURE 5. FACE PIATE WITH BOLT AND SLOT ASSEMBLY 
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The carriage consisted of a U-shaped, welded aluminum platform 
fitted with four 3/4 inch linear bearings. The sides of the platform 
were constructed of two inch OD aluminum pipe fitted at each end with 
aluminum bearing blocks designed to accammodate the linear bearings. 

The structural transverse member of the carriage was a 2 x 3 inch alum- 
num box which was welded to the pipes at one end. The cantilever beam 
used to measure forces on the cylinder was mounted on this transverse 
member. At the opposite end of the aluminum pipe, the bearing blocks 
were connected by a 1/2 inch diameter transverse shaft which accommodated 
the cylinder holding and balancing jig. 

The U-shaped platform was mounted through the linear bearings 
onto two case-hardened 3/4 inch shafts. The shafts were supported and 
held parallel by two six inch channel end pieces which transversed the 
wave channel. The six inch channel end pieces were clamped to the top of 
the wave channel at a position over the observation window. The driving 
mechanism oscillated the carriage sinusoidally. The driving rod was 
connected to the 2 x 3 inch aluminum transverse member by a stiff link 
Which passed through an accommodation hole drilled through the intervening 
Six inch channel end piece. 

3. Cylinder Holding and Balancing Jig 

The cylinder holding and balancing jig consisted of two long arms 
pivoted near their centers as shown in Figure (7). The jig was designed 
to hold the test cylinder at the desired depth of submergence. It was 
designed to have as low a moment of inertia with respect to its pivot 
point as possible in order to maintain a high system natural frequency 


while having adequate sensitivity of force measurement. A counter balance 
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was placed on the amms opposite the cylinder and adjusted so that the 
center of gravity of the system was at the pivot point. As a result of 
this balance adjustment the only force (or moment) acting on the arms 
was that due to the water; all other inertial loads were balanced out. 

A twelve inch long arm was welded at the longitudinal center of 
the center drum and functioned to convert a torque about the rotational 
center of the drum to a linear a This force acted, through a ball 
bearing linkage, on a small cantilever beam. ‘Two attachment points on 
the beam were provided so that a large range of forces could be measured 
without overstressing the beam. The stiffness of the beam was adjusted 
for adequate sensitivity while providing the maximm possible system 
natural frequency. The design natural frequency was approximately 5.0 
times the maximum forcing frequency. 

4, Beam and Link 

The beam and link assembly as shown in Figure (8) functioned to 
rigidly connect the cylinder holding and balancing jig to the carriage. 
Strain gauges were applied to the beam to provide a means of measuring 
forces imposed upon the cylinder and transmitted mechanically through 
the assembly. The strain gauge output was calibrated using dead weight 
loads applied horizontally to the cylinder through a pully arrangement. 

5. Test Cylinders 

Two cylinders as shown in Figure (9) were constructed to provide 
maximum flexibility in testing over a wide range of dimensionless fre- 
quency parameters. A four inch OD cylinder was constructed of thin- 
walled aluminum tubing and a six inch OD cylinder was constructed of 


plexiglass. The cylinders were slotted to accammodate one end of the 
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holding and balancing arms. Bulkheads were placed adjacent to the slots 
to provide an attachment point for the ends of the arms. Silicone sealant 
was used to seal the cylinder against leakage after the arms were attached. 
This allowed the cylinder to be unflooded resulting in a minimum fixed 
mass. 

End pieces as shown in Figure (10) were constructed of plexiglass 
for each cylinder. The end pieces were fitted into the ends of the 
cylinders with 1/8 inch "0" rings to provide the necessary seal. End 
pieces were "dished" on the inner surface to reduce weight. The outer 
surfaces of the end pieces had a modified labyrinth seal cut into the face 
to prevent water from flowing between the wave channel wall and the 
cylinder end when the cylinder was oscillated. 

6. Wave Height Gauge 

The wave height gauge, as shown in Figure (11), consisted of a 
small diameter aluminum rod which acted as the basic frame to which a 
plexiglass mounting table was attached. The elements of a Wheatstone 
bridge were mounted on this table, and two 30 gauge copper wires spaced 
1/2 inch apart were extended from the table and attached to an insulator 
mounted at the end of the aluminum rod. The other end of the rod was 
fitted to a block and stand which allowed the gauge to travel +6 inches 
in the vertical direction. The block and stand provided a convenient, 
accurate means of calibration through adjustment of the submergence of 
the wires. The complete unit was positioned over the channel and 
centered at about five feet from the carriage. 

The wave height gauge was of the parallel wire resistance type. 
The gauge operated on the principle that the conductance between two 


parallel submerged wires varies proportionately to the length of wire 
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submerged. The wave height gauge was connected in parallel to one leg 
of a Wheatstone bridge as shown in Figure (12). To campensate for the 
inability of the carrier preamplifier to offset the resistance of the 
immersed parallel wires, an external fixed resistance was connected 
across an opposite leg of the bridge. The Wheatstone bridge was connected 
with shielded cable to a carrier preamplifier and recorder. 
7. ™&R. P.M. Switch 
The rpm switch as shown in Figure (13) was designed to measure 
the frequency of the motion imposed on the cylinder. At the end of each 
stroke, the switch was closed and the event recorded on the chart. This 
mark was also used to determine the phase shift between the force and 
acceleration. 
8. Carrier Preamplifier and Recorder 
A two-channel Hewlett Packard 8805B carrier preamplifier and 
recorder was used. One channel recorded force on the cylinder and the 
second recorded wave height. The event marker was utilized through an 
external jack to record the frequency of the imposed motion and to 


establish the motion/force phase relationship. 


B. TEST PROCEDURE 
1. System Selection and Optimization 
The first step in the procedure was to select one of two sets of 
holding and balancing arms to provide the desired d/a ratio and to mount 
them on the center drum. Wooden shims were available in 1/2 inch thick- 
ness which allowed placement of either cylinder at two or three radii 
from the wave channel bottan. Once a set of holding and balancing arms 


was chosen, the dimensionless parameter d/a varied with water depth. 
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The next step was to decide whether it was desired to examine the 
lower or higher range of dimensionless frequency parameters. The larger 
six inch OD plexiglass cylinder was used for the lower range of frequency 
parameters because its larger displaced volume produced forces which were 
more readily measurable. However, at higher frequencies the larger 
cylinder resulted in natural to forcing frequency ratios too close to 1.0 
to give valid results. The larger six inch OD cylinder was, therefore, 
used at the lower motor speed range of from 22 to about 100 rpm. The 
smaller four inch OD aluminum cylinder was used at upper speed range of 
from about 100 to 180 rpm. 

Once the combination of holding and balancing arms and cylinder 
was chosen, the cylinder was attached to the arms. The cylinder was 
sealed and the cylinder holding and balancing jig was balanced using the 
adjustable counter balancing weight. After balancing, the jig was placed 
in position on the carriage and the carriage properly positioned and 
clamped to the top of the wave channel. © 

The desired stroke was adjusted at the variable speed trans- 
mission face plate. The stroke was maintained at a value of 0.1 a to 
insure that separated flow did not occur. 

The channel was filled and data runs conducted at water levels 
which were 2, 3, and 4 cylinder radii above the cylinder center. 

After the wave channel was filled to the desired level, the 
Wheatstone bridges were balanced and a test run was made through the 
frequency range to determine the optimum attenuation settings, chart 
speed, and beam position to obtain the most readable results. At this 
time the entire assembly was "plucked" to determine the natural to forced 
frequency ratio. If the value of this ratio was less than five, the beam 


was adjusted to the position which provided the stiffest beam. 
34 





2. System Calibration 


a. Force Calibration 

After balancing the bridge and selecting the appropriate 
attenuator setting, a piece of light weight plastic line was tied loosely 
around the center of the test cylinder and bent around a low friction 
pulley located at the top of the channel, 100 inches from the test cylirder. 
Weights were hung from the line and the resulting deflection on the 
recording paper recorded. The angle which the string made with the hori- 
zontal waS accurately measured and the resultant horizontal component of 
force determined. The strain gauge bridge was calibrated in this manner 
in both directions separately. The results were plotted accurately to 
insure linearity and to be used to interpret data. The calibration of 
the strain gauge bridge was checked frequently between data runs. See 
Appendix A for a typical force calibration plot. 

b. Wave Height Calibration 

After balancing the bridge and selecting the proper attenu- 
ator setting the wave height gauge was calibrated by utilizing the block 
Stand to immerse and withdraw the wire a known distance. The following 
procedure was used. The gauge was raised 1/2 inch simulating a decrease 
in water level, then lowered 1/2 inch below the null, simulating an 
increase in water level. Then, the gauge was raised one inch and lowered 
one inch, etc. This procedure closely simulated the actual passage of a 
wave profile with crests and troughs. Resulting deflections on the 
recording paper were accordingly recorded and the results accurately 
plotted to insure linearity and to be later used to interpret data. The 
calibration of the wave height meter was checked frequently between data 


runs. See Appendix A for a typical wave height calibration plot. 
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c. Event Marker Calibration 
The ability of the event marker to record the frequency of 
the imposed motion and the phase angle was dependent upon the accuracy 
of the chart speed. The accuracy of the chart speed was checked against 
the internal timer of the recorder. No discrepancy was found. 
3. Datamkun 

After system selection and optimization was accomplished and system 
calibration completed, a data run was conducted. The system is shown in 
Figure (14). The data run was commenced with the variable speed trans- 
mission set at the lower value of rpm corresponding to the lower value of 
the range of desired frequency parameters. The inductance motor and 
recorder were turned on simultaneously. The recorder was allowed to run 
until the channel recording wave height indicated a steady state had been 
reached. The recorder was then secured and the transmission adjusted to 
the next rpm value corresponding to the next desired frequency parameter. 
The meter was then sevured and the wave action in the channel allowed to 
subside before the next data point was recorded. This procedure was 
repeated until the desired range of frequency parameters had been covered. 

At the end of each data run, a calibration check was conducted on 


both the wave height gauge and strain gauge bridge. 
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IV. EXPERIMENTAL RESULTS 


As shown in Equation (3) in THEORETICAL CONSIDERATIONS, the force 
coefficient can be represented as a function of relative water depth, 
relative depth of submergence, relative amplitude, and a dimensionless 
frequency parameter. The dynamic forces on the object may be expressed 
as the sum of two components, one component in phase with acceleration 
and one component in phase with the velocity of the cylinder. The 
coefficient of the acceleration component represents added mass and the 
coefficient of the velocity component represents damping. It follows 
that these two coefficients are determinable by measuring the force on 
the cylinder and the phase angle between force and acceleration. 

Additionally, if a control volume were drawn around the system, it 
could be shown through a consideration of conservation of energy that 
the wave energy flux out must equal the energy associated with damping 
if viscous effects are disregarded. Therefore, wave height is directly 
related to the degree of damping. 

Accordingly, a program was designed to measure forces, phase angles 
between force and acceleration, and wave heights within a range of inde- 
pendent dimensionless variable limited by time and equipment available. 

Twelve data sets were taken at six variations of relative water depth 
(d) and relative depth of submergence (h). The six variations were: 
ee 2.0, h = 4.037 cd = 2.0, h = 5.0; d= 3.0, h= 5.0; d = 3.0, h = 6.0? 
d= 4.0, h= 6.0; d= 4.0, h= 7.0. All data sets were taken ata 
relative displacement (x) of 0.1. The range of dimensionless frequency 


parameters (v) for each run was from 0.04 to 1.7. The low frequency half 
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of each run was conducted with a six inch OD plexiglass cylinder and the 
high frequency half with a four inch OD aluminum cylinder. An overlap 
between high frequency and low frequency parameters of approximately 0.36 


was provided near mid-range. 


A. EXAMPLE STRIP CHART READ-OUT 

An example of a strip chart read-out is presented in Figure (15). 
ferethis particular run a = 3, d= 6.0, h= 12.0. The values read fran 
the chart for this data point were wave height (WH) = 0.088 inches, 
force = 0.3525 pounds, period = 143 millimeters, phase angle (PA) = 8 
millimeters, and chart speed = 100 millimeters per second. 

The wave profiles are nearly sinusoidal. The natural frequency of the 
entire mechanism is seen superimposed upon the force profile. It is 
apparent that the natural frequency is considerably greater than the 
forcing frequency. Difficulty was encountered in reading phase angle 
from the chart. A chart speed of 100 millimeters per second was required 
to resolve phase angle. Because of the superposition of natural fre- 
quency upon the force profile, it was necessary to "fair" the force 
profile to properly read phase angle. The result of reading phase angle 
in this manner led to inaccuracies, especially at lower frequencies where 


forces were small and resulting force profiles relatively flat. 


B. DATA REDUCTION 
A computer program was written to reduce experimental data. The 
program is presented in Appendix B with tabulated results for each data 
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C. PRESENTATION OF RESULTS 

Experimental results are presented in graphical form in Figures (16-23). 
For each combination of the dimensionless parameters d and h, a graph of 
force coefficient, phase angle, and wave height ratio was drawn to exhibit 
dependence upon the frequency parameter. In a few cases theoretical curves, 
provided by Garrison in an unpublished paper, were drawn in for comparison. 
Gmooth curves were "faired" in cases where Garrison's theoretical data was 
not camnplete. 

1. Force Coefficient 

Despite the scatter in the values of force coefficient, results 
compared favorably with Garrison's theoretical curves, except at the lower 
end of the frequency spectrum. See Figures (18), (22), and (23). 
A consideration of the boundary value problem associated with the 

problem resulted in the conclusion that the free surface boundary would 
act as a solid boundary at very low frequencies and as a constant pressure 
boundary at very high frequencies. It follows that force coefficient as 
a function of frequency should approach a constant value asymtotically at 
both the high frequency and low frequency end of the data range and vary 
in a continuous manner between. Also, since the constant pressure boundary 
at the high frequency end tends to eliminate pressure differences in the 
fluid, it follows that the constant asymtotic value at the higher frequency 
end should be lower. Further, it follows that as the relative proximity 
of the free surface increases, the effect should be more pronounced. The 
results as presented in Figures (16, 18, 20-23) indicate that these effects 


were observed. 
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Wave Height Ratio 





FIGURE 17. WAVE HEIGHT RATIO 
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FIGURE 19. WAVE HEIGHT RATIO 
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Figure (16) indicates that for d = 2.0 and h = 4.0 the lower 
frequency value of force coefficient was approximately 1.5. The range 
of frequency parameters did not extend to high enough values to deter- 
mine the high frequency asymptote. Theoretical data was not available 
for this figure, therefore a mean line was drawn through the data. Force 
coefficient appeared to reach a minimum value of 0.8 at a frequency 
parameter of 1.15. 

At the greatest depth of submergence, (d = 4.0, h = 6.0), Figure 
(22) indicates that the lower frequency value of force coefficient 
decreased to 1.38. Again, the high frequency asymptotic behavior was 
undeterminable. Garrison's theoretical data at these parameters indicate 
that the force coefficient is still lower than experimentally indicated 
throughout almost the entire range of frequency parameters by as much as 
10%. The minimum value of force coefficient at this larger submergence 
depth was approximately 1.0 at a frequency paremeter of 1.10. 

The effect of the proximity of the impermeable bottcam surface 
can be observed by camparing Figures (16) and (18). During these two 
data runs, relative submergence depth (d) was held constant at d = 2.0 
While relative water depth (h) was varied. Similar comparisons were 
made at d = 3.0 and d = 4.0. In each case, the proximity of the impermeable 
bottam surface appeared to increase the force coefficient over the entire 
range tested. 

It may be noted that the scatter in values of force coefficient 
was greatest at the low frequency end of the range of frequency coeffici- 
ents. This was attributable to the nature of the error in force measure- 
ment. It was cbserved that although the error in force measurement could 


be assumed relatively constant, the magnitude of force measured at values 
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of v< 0.1 was very small, often smaller by an order of magnitude than 
the magnitude of the force measured at values of u> 0.1. It follows 
that with constant absolute error, the relative error can be expected to 
be an order of magnitude larger. In every case scatter was reduced at 
higher frequencies where forces were larger. Also in the area where 
overlap occurred between high and low frequency data sets, the results 
appear consistant and repeatable. It was observed that for values of 

0.1, almost all of the data was consistant within 10%. 

2. Phase Angle and Wave Height 

Values of phase angle and wave height compared favorably with 
Garrison's theoretical curves. See Figures (18), (19), (22), and (23). 
As stated earlier, the added mass coefficient is determinable if the force 
coefficient and the phase angle are known. Further, a theoretical con- 
sideration resulted in the conclusion that added mass coefficient may be 
comouted by multiplying force coefficient by the cosine of the phase angle. 
The balance of the force coefficient is attributable to damping and con- 
sequently, the damping coefficient may be obtained by multiplying the 
force coefficient by the sine of the phase angle. Also, as stated earlier, 
from a consideration of conservation of energy, wave height may be used 
aS a measure of the degree of damping. 

The value of phase angle and wave height at the extremes of the 
frequency spectrum was predicted by intuitively assuming that at very low 
frequencies the force and motion would be in phase, and at very high 
frequencies damping would be non-existant. It was predicted that phase 
angle and wave height would be non-existant where force coefficient 
approached a constant value asymptotically. It followed that force 


coefficient and added mass coefficient would be syncnymous near the 
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asymptotes. The effect of a free surface was predicted to be more pre- 
daninant at lower values of relative depth of submergence. The results 
as presented in Figures (16-23) indicate that these effects were 
observed. 

Figure (16) indicates that for d = 2.0 and h = 4.0 the phase angle 
was zero for frequency parameters greater than 1.5. The lower frequency 
data was not complete for phase. However, Figure (17) indicates that wave 
height was zero for frequency parameters less than 0.4. Figures (16) and 
(17) show phase and wave height reached maximum values of 0.48 radians 
and 0.74 inches respectively. 

At the greatest depth of submergence, Figure (22) indicates that 
for d = 4.0 and h = 6.0 the phase angle was zero for frequency parameters 
greater than 0.9. Again, the phase angle did not vanish at lower fre- 
quencies but wave height vanished for frequency parameters less than 0.4. 
Maximum values were about 0.15 radians and 0.24 inches for phase angle 
and waveheight respectively. As predicted, values were considerably less 
than those observed at a smaller relative submerged depth. 

The effect of the proximity of the impermeable bottom on the phase 
angle can be observed by comparing Figures (16) and (18). The maximm 
value appears to be reduced and shifted slightly to the left. A similar 
reduction and shift to the left were observed in the maximum value of 
waveheight between Figures (17) and (19). 

Again, despite the scatter in the values of phase angle, the 
results memes favorably with Garrison's theoretical curves. In cases 
where theoretical curves were not available the trend in data was "faired 
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All the remarks above pertaining to measuring small magnitude 
forces with constant absolute error apply also to the measurement of 
phase. Accurate measurement of phase angle was especially difficult in 
the lower frequency ranges where force profiles were flat. The exact 
point where the force changed direction was difficult to determine 
accurately. Consequently, more scatter was observed in the phase data. 
It was for this reason that no attempt was made to plot the added mass 
dependency upon frequency parameter. Rather, the resulting sets of force 
coefficient data and phase data were considered separately. Added mass 


is presented in tabular form as computer program output in Appendix B. 
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V. CONCLUSIONS 


An experimental study of the added mass of a right circular cylinder 
under a free surface has been conducted. From the results of this study, 
the following conclusions were drawn. 

1. Results agreed within 10% with theoretical results proposed by 
Garrison in an unpublished paper. 

2. The relative proximity of a free surface induced a pronounced 
variation in added mass coefficient from the classical value of 1.0. 

3. When the cylinder was oscillated at very low frequencies and at 
at very high frequencies, the value of the added mass coefficient was 
constant. This constant value was larger at a low frequency than at 
high frequency due to the fact that at low frequency the free surface 
acts aS a rigid boundary. 

4. The proximity of an impermeable bottcam surface had the effect 
on the value of the added mass coefficient of raising its magnitude at 
all frequencies. 

5. A phase angle between induced force on the cylinder and the 
acceleration experienced by the cylinder was observed to be a function 
of frequency. ‘This phase angle vanished at high and at low frequencies. 

6. The proximity of either a free surface or an impermeable bottom 
Surface had an effect of raising the magnitude of the maximm value of 
phase angle and shifting the occurrence of this maximum to a lower 


frequency value. 
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7. When the cylinder was oscillated, induced waves were observed. 
Wave height was a function of frequency of oscillation. Induced waves 
vanished at low and at high frequencies. 

8. The proximity of either a free surface or an impermeable bottom 
surface had the effect of raising the magnitude of the maximm valve of 
wave height and shifting the occurrence of this maximm to a lower 


frequency value. 
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APPENDIX B COMPUTER PRORAM 


THIS PROGAM WAS USED TO REDUCE EXPERIMENTAL DATA. 
DEFINITION OF PARAMETERS USED ARE AS FOLLOWS: 


T=PERIOD OF FORCING FUNCTION( SECONDS) 


SIGMA=ANGULAR FREQUENCY OF FORCING FUNCTION 
(RADIANS/SECOND) 


G=GRAVITATIONAL CONSTANT CFEET/SECOND**2 ) 
HB=WATER DEPTH(C INCHES) 

AB=CYLINDER RADIUS( INCHES) 

RHO=WATER DENSTITY(SLUGS/FEET**3) 
LB=CYLINDER LENGTH( INCHES) 

DB=DEPTH OF SUBMERGENCE(C INCHES) 
WH=TROUGH TO CREST WAVE HEIGHT (INCHES) 
F=MAXIMUM FORCE(C POUNDS ) 

XB=HALF AMPLITUDE OF MOTIONCINCHES) 
TMM=PERIOD READ FROM CHART(MILLIMETERS) 
PMM=PHASE SHIFT READ FROM CHART( MILLIMETERS) 
CSPEED=CHART SPEED(MILLIMETERS/SECOND) 
PHASE=PHASE SHIFT(RADIANS) 

NCARD=NUMBER OF DATA CARDS 
ANU=DIMENSTONLESS FREQUENCY PARAMETER 
WHR=DIMENSTONLESS WAVE HEIGHT RATIO 

mo OMENS ONESESS FE@RCE CGEFFICIENT 
AM=ADDED MASS COEFFICIENT 

DN=DAMPING COEFFICIENT 
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APPENDIX B (CONTINUED) 


an 
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& 
Me 


"DB=",F8.41 5X_,'HB=" »F8.4, 
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PosmeGhe SPECIFIC VALUES TOG CONSTANTS TO BE LATER USED 


G=32.174 
RHO=62.4/G 


THIS PROGRAM WAS DESIGNED TQ HANDLE TWELVE DATA SETS 
DO 300 K=1,12 

READ IN PARAMETERS WHICH CHARACTERIZE DATA SET 
READ(5,10) AB+sDB,yHB;XByNCARD 

oie ranAMElERS FOR LATER REFERENCE 

WRITE( 6,15) AB,DB,»HB,XByNCARD 

COMPUTE WALWES OF DIMENSIONLESS PARAMETERS 
D=DB/ AB 

K=XB/A6 

PRINT DIMENSTONLESS PARAMETERS 

WRITE(6,20) X;D,H | 

PRINT COLUMN HEADINGS FOR INPUT / QUTPUT PARAMETERS 
WRITE(6,65) 

COMMENCE DATA REDUCTION 

DO 200 T=1,NCARD 

Peap }NeUT “PARAMETERS 

READ(5,25) FyWH, TMM,PMM,CSPEED 

C@mPUTE PERT@D OF FORCING FUNCTION 
T=TMM/CSPEED 

COMPUTE PHASE ANGLE 

PHASE=2.%3.14159%PMM/T MM 

COMPUTE ANGULAR FREQUENCY OF FORCING FUNCTION 
SIGMA=2.*3.14159/T 

COMPUTE DIMENSIGNLESS FREQUENCY PARAMETER 
ANU=SIGMA**2*AB/(G*12) 
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APPENDIX B (CONTINUED) 


COMPUTE WAVE HEIGHT RATIO 
WHR=(WH/2.)/XB 

COMPUTE VOLUME OF CYLINDER 
VOL=3.14159*AB *¥*2%1.25/144. 
Cemeure FOSCE COEFFICIENT 
FC=F/(RHO*VOL*SIGMA*®*2*XB/12.) 
ComreUrTe ADDED MASS CUEFFICTENT 
AM=FC*COS (PHASE) 

CGMPUTE DAMPING COEFFICIENT 
DN=FC*S INC PHASE) 

CONRUTATLON COMELE EG FRTAT, VAs 
AT THE NEXT DATA POINT IN THIS 


WRITE(6,;80) SIGMA, F,CSPEED,ANU,FC,PHASE,;AM,DNs WHR 
CONTINUE 


CORO NEv DATA SET. 
CONTINUE 

Soe 

END 


UETORTPAmRAMETERS AT 
Se Om eeePe 10 EOSK 
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APPENDIX B (CONTINUED) 
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APPENDIX 8B (CONTINUED) 


DATA RUN NUMBER 2 
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APPENDIX 8 (CONTINUED) 


DATA RUN NUMBER 3 
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APPENDIX B (CONTINUED) 


DATA RUN NUMBER 4 
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